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A maximum likelihood constrained time-domain adaptiv'e beam- 
former was evaluated for its detection performance using short-period data 
from the Korean Seismic Research Station (KSRS), A total of 132 events in 
Eurasia was used for the detection study, and two hour-long noise samples 
were used for false alarm investigation. The results indicated that the pre- 
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SECTION I 


INTRODUCTION 

A. THE ADAPTIVE BEAMFORMING SYSTEM 

The adaptive beamformine (ABF) system of Barnard (1974) is 
a maximum likelihood multichannel time- domain adaptive filtering beamform- 
er. The design goal is to minimize the filter output power subject to unity 
response constraints in the desired steer direction. The constraints are de- 
signed for passing a plane wave signal arriving from the desired steer direc- 
tion in ar ^informing. 

■ e ABF output y(t) is formed by applying a convolution filter 
to eac;. cj^annel and summing the outputs of all channels: 

k’. N 

yft) = ^ ^ a.(j) X.(t-j) 

i=l j = -N 

th 

v^'here a^(j) is the filter weight for the i channel at a lag of j sample 
points, x^(t-j) is the value of the channel i at time t-j. k! is the number of 
channels, and 2N + 1 is the total length of the filter in sample points. Prior 
to forming the filter outputs each channel is time- shifted to time-align energy- 
arriving from the desired steer direction. 

The adaptive filter weights a^(j) are updated by the following 

algorithm : 

a"^'^(J) = a"^'*(j) ^ A(t)y(t) j^x(t-j) - 


I-l 


L. 


J 


where 




and X(t) , which controls the adaptation rate at time t , is computed as fol- 
lows: 


X(t) = 


M N 


i=l j = -N 


For this algorithm, H is the convergence rate discussed in the report and is 
in general in the range 0 < 4 < 2. The constraints imposed on the filter 


weights a.(j) are: 


= 1 


for j ^ output point 


for j = output point. 


OBJECTIVES 


The objectives of this study are as follows: 

• To determine the adaptive noise reduction relative to beam- 
steering and to determine the false alarm probability for the 
adaptive beamforming detector using short- period data from 
KSRS. 

• To determine the ABF detector's detection performance and to 
estimate the increased array detection capability by use of the 
ABF detector. 

Section II presents a brief description of the short-period ar- 
ray at the Korean Seismic Research Station and the data base used in this 
study. Section III presents the results in terms of the false alarm study, 
detection probability measurements, simulation study, and short- period ar- 
ray response. Maj^r results 'and'conclusions are given in Section IV. 
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SECTION II 
DATA BASE 


A. ARRAY 

The short-period array at the Korean Seismic Research Station 
(KSRS) is a 19-element hexagonal array with vertical components. It has an 
aperture of about 10 km with two rings (for configuration, see Prahl et al. , 
1975). 

B. DATA 

Data used in this work were taken from the periods of Novem- 
ber 1974 and January and February 1976. Two hour-long noise samples were 
used for noise reduction and false alarm study. A total of 132 events, 38 of 
them in central Eurasia (November 1974) and 94 of them in Kurile- Kamchatka 
(January and February 1976), were used for the detection study. 






r 

SECTION m 
RESULTS 



Subsection A presents noise study and false alarm probabilitv 
measurements. Subsection B presents the results of detection probability 
when the ABF was applied to signals. In Subsection C, a simulation study 
was conducted by burying a scaled signal in noise, followed by array responses 
in Subsection D. 

A. NOISE STUDY AND FALSE ALARM PROBABILITY 

Two hour-long noise samples in November 1974 were process- 
ed to study the ABF noise reduction relative to beamsteering. In order to find 
the optimum performance for the ABF, a number of prefilters and two differ- 
ent ABF filter lengths (15 and 31 points) were applied to the single channel 
data. The prefilters applied were wideband (unfiltered), low frequency pass- 
band filter (Figure III-l), array beam optimum detection filter (Figure III-Z), 
single sensor optimum detection filter (Figure III- 3), and noise whitening fil- 
ter (Figure I1I-4). The filter in Figure III-l was implemented in the ABF sys- 
tem (originally intended for long-period data), and ABF v.dll default to this fil- 
ter. The filter in Figure III-2 was designed on the basis of output beam signal 
and noise spectra to optimize the signal -to- noise ratio for the array beam and 
therefore is called array beam optimum detection filter. For the same purpose, 
Figure III- 3 was designed from the averaged single- sensor and noise spectra. 
Finally, the filter in Figure 1II-4 was designed from noise spectra only to 
whiten the noisi . 

Table III-l shows the processing results from two noise samples 
(days 312 and 325, 1974), The array steer direction was zero degrees wdth 15 
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FIGURE III- 2 

FREQUENCY RESPONSE FOR REAM OPTIMUM 
DETECTION FILTER 
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ADAPTIVE DEAMFORMER NOISE REDUCTION RELATIVE TO BE AMSTEE RING 



Numbers in this table arc 20 Log (RMS noise of beamsteer/RMS noise of ABF) decibels. 












km per second velocity. The processing gains in the table suggest that the 
ARF achieved the best noise reduction with the low frequency passband filter, 
because the low^ frequency noise is mostly surface wave modes which are 
better correlated than the noise in higher frequency ranges which are random 
and compres sional modes (Prahl, et al. , 1975). Noise spectral analysis 
(Figure III-5) suggests that the ABF performance decreases as frequency in- 
creases. 

Noise outputs from the beamsteer and the ABF beams W'ere 
used for conventional power detector computation. The ratio of the short- 
term average to the long-term average of these outputs in decibels was form- 
ed, where a 30 second sliding window was used for long-term average com- 
putation, while various time gates of 0. 8, 1. 6, 3. 2, and 6.4 seconds were 
used for short-term' averages. The mean and standard deviation were com- 
puted from these output ratios. Table III- 2 and Table III- 3 summarize the 
results of these calculations. It is noteworthy from these results that there 
is very little difference between beamsteer and ABF results. Also, there is 
little day-to-day change in the results obtained for the winter noise. 

B. DETECTION PROBABILITY MEASUREMENTS 

In the following measurements, predicted event azimuth and 
velocity for beamsteering and ABF were computed from the event locations 
listed in the Norwegian Seismic Array (NORSAR) biilletin. Signal-to-noise 
ratios (SNR's) were computed using events with high SNR on single-channel 
data. The signal-to-noise ratio was the peak-to-peak signal amplitude 
divided by two- minute RMS noise amplitudes in decibels taken prior to signal 
arrival. For strong signals, the ABF SNR impro\ement relative to beam- 
steering was about to 8 dP using the low frequency passtsand filter. The 
improvement uas due to noise reduction achiev'ed l>y the AI^F consistent witli 
Tatile III- I. When the low frequency passband prefilter was applied to the 
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TABLE III -2 


MEAN AND STANDARD DEVIATION FOR ABF AND BEAMSTEER 
DETECTORS (NOISE SAMPLE: DAY 312 HOUR 15:00 TO 15:57, 
1974; 15-POINT ABF LENGTH; 0,5 CONVERGENCE RATE) 



Single- 

Channel 

Prefilter 


Unfiltered 


Single-Sensor 

Optimum 

Detection 

Filter 


Low 

F requency 
Band Filter 
(0. 5-1.1 Hz) 


Array-Beam 

Optimum 

Detection 

Filter 


Whitening 

Filter 


Gate 

(seconds) 


Beamsteer Output (dB) 

ABF Output (dB) 

Average 

Standard 
Deviation (o) 

Average 

Standard 
Deviation (a) 




0. 8 

- 

1. 6 

- 

3. 2 

- 

6.4 

- 

0. 8 

- 

1 . 6 

- 

3, 2 

- 

6.4 


0, 8 


1.6 


! 3, 2 

1 


' 6. 4 

- 


0.435 j 
0.218 

2.003 

1.422 

— 



2. 073 
1. 447 
0. 838 
0.414 



-1,464 
- 1 . 927 
-0. 552 
7. 308 



4. 770 
3. 91 1 
2. 917 
2. 014 


3. 777 
2. 940 
2. 252 
1. 699 


. 246 
0. 723 


■ 1 . 944 
■1. 239 
■0. 685 
•0. 331 


•1. 248 
■0. 73P 
•0.420 


4. bOO 


2, t.54 



3. 500 
2. b35 
1 . 9b6 































TABLE III- 3 




MEAN AND STANDARD DEVIATION FOR ABF AND BEAMSTEER 
DETECTORS (NOISE SAMPLE: DAY 325 HOUR 05:05 TO 05:57 
1974; 15- POINT ABF LENGTH: 0.5 CONVERGENCE RATE) 


Single- 


Beamsteer Output (dB) 

ABF Output (dB) 

Channel 

Prefilter 

(seconds) 

Average 

M 

Standard 
Deviation (ct) 

Average 

(A) 

Standard 
Deviation (a) 


0. 8 

-2. 196 

4. 775 

-1.492 

3.830 

Unfiltered 

1.6 

-1. 313 

3. 690 

-0. 875 

2. 926 


3. 2 

-0. 654 

2. 529 

-0. 506 

2. 210 


6, 4 

-0. 326 

1. 750 

-0. 275 

1. 628 

Single-Sensor 

0. 8 

-1. 332 

3. 630 

-1. 295 

3. 593 

Optimum 

Detection 

1. 6 

-0. 800 

2. 764 

-0. 763 

2. 710 

Filter 

3. 2 

-0. 442 

2. 053 

-0. 418 

1. 977 


6. 4 

-0. Z32 

1.496 

-0. 215 

1.428 

Low 

F requency 

0.8 

-2, 042 

4. 775 

-1. 900 

4. 549 

1. 6 

-1.433 

3. 923 

-1. 228 

3. 593 

Band I ilter 
(0. 5- 1 . 1 Hz) 

3. 2 

-0. 856 

2. 995 

-0. 700 

2.697 


6.4 

-0.431 

2. 068 

-0. 346 

1. 936 

Array - Beam 

0. 8 

-1. 373 

3. 689 

-1.110 

3. 322 

Optimum 

Detection 

1. 6 

-0. 817 

2. 803 

-0. 624 

2.444 

F i Iter 

3. 2 

-0. 432 

2. U3b 

-0. 338 

1.796 


6, 4 

-0. 213 

1. 438 

-0. 169 

1. 285 


0. 8 

- 1. 532 

3. 875 

-1 . 258 

3. 552 

Whitening 

1. 6 

-0. 991 

3. 077 

-0. 742 

2. 700 

F i Ite r 

3. 2 

-0. 609 

2. 420 

-0. 408 

1. 998 


6. 4 

-0. 348 

1. 863 

-0. 216 

1.478 
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low amplitude signals, the ABF SNR gains degraded to about the beamsteer 
detection performance. 

To estimate roughly the detection capability for the ABF de- 
tector, an ensenible of 38 events in the central Eurasian region in November 
1974 was processed. Single channel data were examined for all events. 
Beamforming was performed for those events for which the 'expected' signals 
were not visible in the single channel data. The criteria for detection was 
that the peak signal amplitude was at least 1. 5 times higher than the peak 
noise in two minute gates preceding the signal arrival time. Gaussian pa- 
rameters for detection probability in terms of bodywave magnitude are given 
in Table 111-4. In general, the ABF detector yielded the same detectability 
as the beamsteer detector. 

After receiving the new data from the KSRS array, a total of 
94 events in the Kurile -Kamchatka region were processed using the January 
and February 1976 data. Two convergence rates (0.05 and 0.5) were used 
with the 15 -point filter length and a detection was claimed for the AFB if 
either or both of the two convergence rate beams yielded a detection. For 
the low frequency passband, the ABF update was stopped based on apparent 
SNR output criteria greater than 3 dB. 

Figure II1-6 presents the least- squares estimate of detection 
probability for the beamsteer detector (Shen, 1976). The upper part of the 
figure shows the histogram of detections and the lower part the Gaussian 
probability curve which indicates a 50 percent detectable bod^nvave magnitude 
of 4. 12 for the beamsteer detector. The lower magnitude capability for 
Kurile- Kamchatka as compared with central Eurasia is due to its smaller 
epicentral distance of 1 7 degrees. The results in the low frequency passband 
and in the single- sensor optimum detection passband for the beamsteer and 
the ABF detectors are presented in Table III-5. 
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TABLE III -5 


I GAUSSIAN PARAMETERS OF DETECTION PROBABILITY 

^ FROM LEAST-SQUARES ESTIMATE (KURILE- 

KAMCHATKA, JANUARY AND FEBRUARY 
, 1976. 94 EVENTS) 


Processor 

Filter 

Mean (/u) 
50 Percent 
Detectable 

(mb) 

i 

1 

Standard j 

Deviation (a) 1 

'"b* i 

ABF 

High 

F requency 
Passband='= 

4. 10 

i 

0. 35 ; 

1 

1 

Low 

F requency 
Passband 

4, 42 

1 

0. 68 ! 

Beamsteer 



High 

F requency 
Passband 

4, 12 

0. 34 

Low 

F requency 
Passband 

4, 54 

0. 63 


Sintzlr- sensor optimum detection filter. 
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By stopping the ABF update upon the signal-to-noise ratio 
greater than 3.0 dE, signal degradation decreased somewhat, but false 
alarms increased considerably. This was due to the fact that with the high 
convergence rate used in the processing the ABF output amplitude increased 
immediately as the ABF update was stopped. For the beamsteer detector, 
the 50 percent detectability was about 0.4 m^ units higher using low frequency 
passband. For the ABF detector, the difference was about 0. 3 m^ units be- 
tween the two passbands used. The estimates in Table III-5 suggest that 
there would not be too much improvement in detection capability at the KSRS 
short-period array by using ABF as compared to beamsteering. 

C. SIMULATION STUDY 

A simulation study was conducted by burying a scaled signal in 
noise. This was done by summing two samples, one containing a signal and 
the other being noise only, in single channel data to form a composite sample 
from which the ABF was designed. The signal-to-noise ratios for the input 
signals were varied in the single- sensor level by scaling the signal sample 
by a factor before adding it to the noise sample. In the simulation, the event 
from Novaya Zemlya was the signal sample. Fifteen sites (sites 1 to 15) and 
a 15-point ABF filter length were used. The signal was buried in the noise 
at various signal-to-noise levels and was processed by both the ABF and 
beamsteer. 

The output signal-to-noise ratios were computed by taking the 
ratio of signal peak-to-peak amplitude to the RMS noise amplitude two minutes 
prior to the signal arrival. Figure III- 7 shows results using single- sensor 
optimum detection filter. The ABF SNR gain relative to beamsteering was 
about 2 dp at hich SNR's. This gain degraded to zero at aV^out 15 dB input, a 
level at which the signal was still easily detectable ijy beamsteering. Below 
15 dP input, this gain was negative. 
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P)eam Output Signal-tcj- Noise Ratio (clB) 



FIGURE III- 7 


SIMULATION RESULTS USLNC 15-POINT FILTER, 0.5 CONVERGENCE 
RATE, AND SINGLE-SENSOR OPTIMUM DETECTION FILTER 

(See Figure III- 3 ) 
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The same simulation was performed using the low frequency 
passband filter. Figure III- 8 presents the results. Two cases were conduct- 
ed: one with continuous update of the ABF and the other with stopping update 
of the ABF when the signal-to-noise ratio is greater than 0 dB. For the form- 
er case, we essentially obtained the same performance pattern as in Figure 
III-7. For the latter case, the ABF SNR gain was always greater than the 
beam steering. 


D. ARRAY RESPONSE 

The detector's response pattern for the KSRS short-period data 
was studied for an event from Greece. The measurements were done by beam- 
forming the event data using various azimuths and velocities. Signal-to-noise 
ratios were computed by taking the ratio of signal peak-to-peak amplitude to 
the RMS noise amplitude two minutes prior to the signal arrivals on the beam- 
steer and the ABF beams, respectively. Figure III- 9 shows the response pat- 
tern for the beamsteer and the ABF detectors. The dots indicate the source 
location in velocity- azimuth space. The ABF signal degradation was more 
pronounced than the beamsteer processor because ABF was not able to take 
the advantage of signal coherence when the single channels were aligned in- 
correctly. 
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SECTION IV 

MAJOR RESULTS AND CONCLUSION 



The adaptive beamforming detector (ABF) was evaluated for its 
performance by use of two hour-long noise samples and a total of 132 events 
in the Eurasian region as recorded at the Korean Seismic Research Station. 
Major results are presented here. 

• For noise reduction, the ABF achieved about a 6 to 10 dB gain 
relative to beamsteering between 0. 5 and 1. 1 Hz and a 2 to 4 dB 
between 1, 5 and 2. 4 Hz. 

• For signal detection, using 38 events from central Eurasia, the 

detection probability estimate yielded a 50 percent detectable 

bodywave magnitude of 4. 35 m for beamsteering and 4. 24 m 

b ^ b 

for the ABF. For 94 events in Kurile- Kamchatka, the 50 per- 
cent detectable bodywave magnitude was 4. 12 m for beam- 

b 

steering and 4. 10 m, for the ABF. It is concluded that the 

b 

ABF provides no improvement over beamsteering for this ap- 
plication. 

• A simulation study suggested that signal degradation by the ABF 
was increased at low input signal-to-noise ratios. Stopping the 
ABF update when signal-to-noise ratio exceeds certain levels 
prev'ented the degradation somewhat at the expense of an in- 
creased false alarm rate. 

The only conclusion which can be drawn from tliese results is 
that the ABF algorithm in its present form as tested here did not yield any 
advantage over a conventional bcamformer for the KSR.S short-period data. 
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